Search for B° Decays to Invisible Final States at Belle 



C. -L. Hsu, 31 P. Chang, 31 I. Adachi, 9 H. Aihara, 48 K. Arinstein, 2 D. M. Asner, 37 V. Aulchenko, 2 T. Aushev, 15 
A. M. Bakich, 43 B. Bhuyan, 10 M. Bischofberger, 29 A. Bondar, 2 G. Bonvicini, 53 A. Bozek, 32 M. Bracko, 24 - 16 

T. E. Browder, 8 M.-C. Chang, 5 Y. Chao, 31 V. Chekelian, 25 A. Chen, 30 P. Chen, 31 B. G. Cheon, 7 K. Chilikin, 15 
I.-S. Cho, 55 K. Cho, 19 Y. Choi, 42 J. Dalseno, 25 * 45 J. Dingfelder, 1 Z. Dolezal, 3 Z. Drasal, 3 D. Dutta, 10 S. Eidelman, 2 
D. Epifanov, 2 S. Esen, 4 H. Farhat, 53 J. E. Fast, 37 V. Gaur, 44 N. Gabyshev, 2 R. Gillard, 53 Y. M. Goh, 7 J. Haba, 9 
T. Hara, 9 K. Hayasaka, 28 H. Hayashii, 29 Y. Horii, 28 Y. Hoshi, 46 W.-S. Hou, 31 Y. B. Hsiung, 31 H. J. Hyun, 21 
T. Iijima, 28 ' 27 K. Inami, 27 A. Ishikawa, 47 R. Itoh, 9 M. Iwabuchi, 55 Y. Iwasaki, 9 T. Iwashita, 29 T. Julius, 26 

J. H. Kang, 55 T. Kawasaki, 34 H. Kichimi, 9 C. Kiesling, 25 H. J. Kim, 21 H. O. Kim, 21 J. B. Kim, 20 
J. H. Kim, 19 K. T. Kim, 20 Y. J. Kim, 19 B. R. Ko, 20 P. Kodys, 3 S. Korpar, 24 ' 16 P. Krokovny, 2 T. Kuhr, 18 
A. Kuzmin, 2 P. Kvasnicka, 3 Y.-J. Kwon, 55 S.-H. Lee, 20 J. Li, 41 Y. Li, 52 J. Libby, 11 Y. Liu, 4 Z. Q. Liu, 12 

D. Liventsev, 15 R. Louvot, 22 K. Miyabayashi, 29 H. Miyata, 34 Y. Miyazaki, 27 G. B. Mohanty, 44 A. Moll, 25 ' 45 
N. Muramatsu, 39 E. Nakano, 36 M. Nakao, 9 Z. Natkaniec, 32 C. Ng, 48 S. Nishida, 9 O. Nitoh, 51 T. Ohshima, 27 

S. Okuno, 17 S. L. Olsen, 41 ' 8 G. Pakhlova, 15 C. W. Park, 42 H. Park, 21 H. K. Park, 21 T. K. Pedlar, 23 R. Pestotnik, 16 
M. Petric, 16 L. E. Piilonen, 52 M. Ritter, 25 M. Rohrken, 18 S. Ryu, 41 H. Sahoo, 8 Y. Sakai, 9 S. Sandilya, 44 
\q \ T. Sanuki, 47 O. Schneider, 22 C. Schwanda, 13 A. J. Schwartz, 4 K. Senyo, 54 M. E. Sevior, 26 M. Shapkin, 14 

C. P. Shen, 27 T.-A. Shibata, 49 J.-G. Shiu, 31 B. Shwartz, 2 A. Sibidanov, 43 F. Simon, 25 ' 45 J. B. Singh, 38 

P. Smerkol, 16 Y.-S. Sohn, 55 E. Solovieva, 15 S. Stanic, 35 M. Staric, 16 M. Sumihama, 6 T. Sumiyoshi, 50 
G. Tatishvili, 37 Y. Teramoto, 36 K. Trabelsi, 9 M. Uchida, 49 T. Uglov, 15 Y. Unno, 7 S. Uno, 9 P. Urquijo, 1 
Y. Usov, 2 S. E. Vahsen, 8 P. Vanhoefer, 25 G. Varner, 8 V. Vorobyev, 2 P. Wang, 12 M. Watanabe, 34 Y. Watanabe, 17 
K. M. Williams, 52 E. Won, 20 H. Yamamoto, 47 Y. Yamashita, 33 Z. P. Zhang, 40 V. Zhilich, 2 and V. Zhulanov 2 

(The Belle Collaboration) 
1 University of Bonn, Bonn 

C^) ' 2 Budker Institute of Nuclear Physics SB RAS and Novosibirsk State University, Novosibirsk 630090 

3 Faculty of Mathematics and Physics, Charles University, Prague 
* University of Cincinnati, Cincinnati, Ohio 45221 
5 Department of Physics, Fu Jen Catholic University, Taipei 
, 6 Gifu University, Gifu 

I • 7 Hanyang University, Seoul 

' University of Hawaii, Honolulu, Hawaii 96822 

9 High Energy Accelerator Research Organization (KEK), Tsukuba 
(N- Indian Institute of Technology Guwahati, Guwahati 

11 Indian Institute of Technology Madras, Madras 
, 12 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 

• <"H ■ 13 Institute of High Energy Physics, Vienna 

1 1Jt Institute of High Energy Physics, Protvino 

, 15 Institute for Theoretical and Experimental Physics, Moscow 

16 J. Stefan Institute, Ljubljana 
11 Kanagawa University, Yokohama 
13 Institut fur Experimentelle Kernphysik, Karlsruher Institut fur Technologie, Karlsruhe 
19 Korea Institute of Science and Technology Information, Daejeon 
20 Korea University, Seoul 
21 Kyungpook National University, Taegu 
22 Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 
23 Luther College, Decorah, Iowa 52101 
2Jl University of Maribor, Maribor 
25 Max- Planck- Institut fur Physik, Miinchen 
26 University of Melbourne, School of Physics, Victoria 3010 
27 Graduate School of Science, Nagoya University, Nagoya 
28 Kobayashi-Maskawa Institute, Nagoya University, Nagoya 
29 Nara Women's University, Nara 
30 National Central University, Chung-li 
31 Department of Physics, National Taiwan University, Taipei 
32 H. Niewodniczanski Institute of Nuclear Physics, Krakow 
33 Nippon Dental University, Niigata 



2 



3J, Niigata University, Niigata 
35 University of Nova Gorica, Nova Gorica 
3G Osaka City University, Osaka 
37 Pacific Northwest National Laboratory, Richland, Washington 99352 
38 Panjab University, Chandigarh 
39 Research Center for Electron Photon Science, Tohoku University, Sendai 
J '° University of Science and Technology of China, Hefei 
J>1 Seoul National University, Seoul 
^ 2 Sungkyunkwan University, Suwon 
^ 3 School of Physics, University of Sydney, Sydney, New South Wales 2006 

^ Tata Institute of Fundamental Research, Mumbai 
45 Excellence Cluster Universe, Technische Universitat Miinchen, Garching 
46 Tohoku Gakum University, Tagajo 
^ Tohoku University, Sendai 
^Department of Physics, University of Tokyo, Tokyo 
^ 'Tokyo Institute of Technology, Tokyo 
50 Tokyo Metropolitan University, Tokyo 
51 Tokyo University of Agriculture and Technology, Tokyo 
r ' 2 CNP, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
53 Wayne State University, Detroit, Michigan 48202 
4 Yamagata University, Yamagata 
5 Yonsei University, Seoul 

We report a search for B° decays into invisible final states using a data sample of 657 x 10 6 
BB pairs collected at the T(4S) resonance with the Belle detector at the KEKB e + e~ collider. 
The signal is identified by fully reconstructing a hadronic decay of the accompanying B meson and 
requiring no other particles in the event. No significant signal is observed, and we obtain an upper 
limit of 1.3 x 10 -4 at the 90% confidence level for the branching fraction of invisible B° decay. 

PACS numbers: 13.20.He,12.15. Ji,12.60. Jv 



In the standard model (SM), the decay B° — > vv pro- 
ceeds through the three annihilation diagrams shown in 
Fig. HJa). This decay is highly helicity suppressed with 
an expected branching fraction at the 10 -20 level [![. 
Because neutrinos participate only in weak interactions, 
the experimental signature is missing energy and momen- 
tum corresponding to the presence of a B° meson in the 
event. New particles hypothesized by physics beyond the 
SM, such as i?-parity violating supersymmetry, can be in- 
volved in these B decays, resulting in a final state with 
only weakly interacting particles and providing the same 
signature as in B° — > vv. For instance, Ref. [2j discusses 
the B decay into a neutrino and a neutralino (x?), shown 
in Fig. [ljb); the branching fraction could be as high as 
10~ 6 — 10~ 7 . Therefore, signals of invisible B decays in 
current B factory data would indicate new physics. So 
far no such signals were observed. The first experimental 
result was provided by the BaBar Collaboration, with 
B(B -» invisible) < 2.2 x 1CT 4 at the 90% confidence 
level [H with a semileptonic tagging method; recently, 
the upper limit was pushed to 2.4 x 10~ 5 with more data 
and improved tagging efficiency by BaBar [4(. 

In this paper we report the result of a search for B de- 
cays to an invisible final state based on the data collected 
with the Belle detector at the KEKB asymmetric-energy 
(3.5 on 8 GeV) e + e~ collider [5J. The data sample con- 
sists of 657 x 10 6 BB pairs accumulated at the T(4S') 
resonance, corresponding to an integrated luminosity of 
606 fb _1 , and an additional 68 fb _1 of off-resonance data 




FIG. I: Feynman diagrams for the SM process via B — > vv 
(a) and for new physics via B° —> Xi~v (b)@|. 



recorded at a center-of-mass (CM) energy about 60 MeV 
below the T(AS) resonance. The Belle detector consists 
of a four-layer silicon vertex detector, a 50-layer central 
drift chamber (CDC), time-of- flight scintillation counters 
(TOF) , an array of aerogel threshold Cherenkov counters 
(ACC), and a CsI(Tl) electromagnetic calorimeter(ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. Outside the coil, the K\ 
and muon detector (KLM), composed of resistive plate 
counters, detects K® mesons and identifies muons. The 
detector is described in detail elsewhere [6j . A GEANT3- 
based . 7] Monte Carlo (MC) simulation of the Belle detec- 
tor is used to optimize the event selection and to estimate 
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the signal efficiency. 

Since the T(4S I ) decays to AA pairs, invisible B de- 
cay candidates are identified by fully reconstructing a 
B meson (A tag ) following the procedure of Ref. Q in 
hadronic modes, and then examining whether there are 
any other particles in the event. The neutral A tag can- 
didates are reconstructed through B° —> D^*^h + de- 
cays, where h + denotes 7r + ,p + ,a^, or Di*^ + [9]. Can- 
didate D* s ^ mesons are identified through the channels 

D*+ -> D+j and D*~ -> if ix~ . Candidate D (s) 
mesons are reconstructed using the following final states: 

K-TT+TT+, A'-TT+TT+TT , K+ K~ 7T+ , K%TT+ , A^TT+TT , and 

- for A+; K-n + , K%tt , K+K~, K-jt+tt , 
and K + K-n + for Df 



and K^tt+tt-it for A ; and 



KgTT + lT ' 7T 
K° S TT+TT-, 

K°K+,K + n-n 

Charged kaons and pions are identified using specific 
ionization from the CDC, time-of- flight information from 
the TOF, and Cherenkov light yield in the ACC. This 
information is combined to form a K-tt likelihood ratio 
T^-K/ir = £k /(£k + where Lk (£„•) is the likelihood 
that the track is a kaon (pion). Tracks with IZk/tt > 0.6 
are regarded as kaons and TZr/tt < 0-4 as pions. The 
typical selection efficiency for a 1.0 GeV/c kaon (pion) 
is 83% (90%) while the misidentification probability for 
1.0 GeV/c kaons (pions) as pions (kaons) is around 6% 
(12%). Neutral Kg — > n + ir~ candidates are identified by 
pairing two opposite-sign charged tracks, both treated as 
pions, and then requiring that this pair have an invariant 
mass near the nominal K§ mass with a vertex displaced 
from the e + e~ interaction point. Candidate A^'s are 
selected from KLM hit patterns that are not associated 
with any charged track [10j . Neutral pions are identified 
using the ir° — > 77 decay and requiring each photon to 
have a minimum energy of 50 MeV and 77 mass between 
0.115 GeV/c 2 and 0.156 GeV/c 2 . The p+ and af meson 



candidates are reconstructed using the ytr 



and 



TT IT TT 



channels. 



The selection of A tag candidates is based on two 
kinematic variables: the beam-energy constrained mass 
M bc = V ^bcam — Pb an d tne energy difference AE = 
Eb — Abeam, where Eb and ps are the reconstructed 
energy and momentum of the A tag candidate in the 
e + e~ CM frame, and Ab ca m is the beam-energy in this 
frame. The B tag candidates are required to have M bc > 
5.22GeV/c 2 and \AE\ < 0.3 GeV. Within this region, 
we define the signal region: 5.27 GeV/c 2 < M\, c < 5.29 
GeV/c 2 and -0.08 GeV < AA < 0.06 GeV. Figure H 
shows the Mbc and AA distributions of the -Btag candi- 
dates in data. If there are multiple A tag candidates in 
an event, the candidate with the smallest \ 2 is retained, 
where x 2 is computed using AA, the D meson mass, and 
the mass difference between the D* and D (for candi- 
dates with a D* in the final state), weighted using their 
expected resolutions. We reconstruct 9.5 x 10 5 neutral 
Atag candidates in total. After identifying the -B tag can- 
didate, we require no additional charged tracks nor 7r° or 
A£ candidates in the rest of the event. 
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FIG. 2: The M bc (a) and AA (b) distributions for the -B tag 
candidates. Candidates having Mbc and AA within the signal 
regions between the two arrows are used to search for B decays 
to invisible final states. 



The dominant backgrounds are from e + e~ — > qq (q — 
u, d, s, c) continuum events and BB decays with a b — > c 
transition (generic B background). Two variables are 
used to distinguish the signal and continuum events: 
cos 9b, defined as the cosine of the angle between the 
Atag flight direction and the beam axis in the CM frame, 
and cos#Tj the cosine of the angle of the A tag thrust axis 
with respect to the beam axis in the CM frame. Clear 
differences in the distribution of each variable between 
signal and continuum background are shown in Fig. [31 
using the MC simulation. We define the fit region as 
-0.9 < cos6>b < 0.9 and -0.6 < cos(9 T < 0.6. The 
variable cos 8b is used in the fit to extract the signal 
yield. Other backgrounds, such as rare A decays via 
b — >• q (q = u,d,s) processes and e + e~ — » t + t~ tran- 
sitions, are also considered in the signal extraction and 
studied using large MC samples. The t + t~ background 
is small and has an event topology similar to the contin- 
uum; therefore, the continuum and t + t~ backgrounds 
are combined and called the non-A background. 

The most powerful variable to identify B decays into 
the invisible final state is the residual energy in the ECL, 
denoted Aecl , which is the sum of the energies of ECL 
clusters that are not associated with the A tag daughters. 
To further suppress the background, minimum energy 
thresholds are required for clusters located in various 
ECL regions: 50 MeV for the barrel (32.2° < 6 < 128.7°), 
100 MeV for the forward endcap (6 < 32.2°), and 150 
MeV for the backward endcap (9 > 128.7°). 

The signal yield for invisible B decays is extracted 
from an extended unbinned maximum likelihood fit to 
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FIG. 3: Normalized distributions of (a) cos 6b and (b) cos St 
for the signal and continuum backgrounds. The solid his- 
togram is the signal and the dashed histogram is the contin- 
uum background. 



the -Eecl and cos 9b distributions. The likelihood is 

- v 
e ^ 3 



N 



c 



II 5I^j(^ECL,COS0 B ) 



(1) 



1=1 



where i is the event identifier; nj is the yield for category 
j, which corresponds to either signal, generic B, rare B or 
non-B background; and ^(Eecl, cos 9b) is the product 
of the probability density functions (PDFs) 'P(Eecl) and 
V (cos 9b), since we have verified that Eecl and cos 9b 
are uncorrelated for each component. For each category, 
the Eecl PDF is modeled as a histogram function, while 
the cos 9b PDF is described by a first or second order 
Legendre polynomial. The non-E? Eecl PDF is con- 
structed from off-resonance data, while all other PDFs 
are obtained using MC simulations. The normalization 
of the rare B background category is estimated from the 
MC simulation and is fixed in the fit. 

The .Eecl simulation is validated using doubly tagged 
events in which the -B tag is fully reconstructed as de- 
scribed above and the other B is identified as B° — > 
D^-e+v (l = e,/i). Candidate D* 



v [l = e,u), Oancudate L) mesons are recon- 
structed via D*~ -> 75°7r~, followed by E>° -> K + u~ 
while D~ is identified as D~ — > K*lir~ and K + tt 



TV 



The track and 7r° selections are applied here. Background 
contributions in the doubly tagged sample are found to 
be negligible; therefore, only loose selections on D and 
D* masses and the mass squared of the undetected parti- 
cles m 2 miss = |Pboam - Ps tag - Pd<»)-£+| 2 (where P de- 
notes the four-momentum of the e + e~ system, Stag, or 
the D^~£ + system) are applied. 

The observed Eecl distributions for doubly tagged 
events, shown in Fig. |H are found to be in good agree- 
ment with MC simulations. The signal yields for control 
modes are obtained by fitting the Eecl spectra while the 



efficiencies are estimated from MC samples. The mea- 
sured branching fractions with their errors, listed in Ta- 
ble U agree well with the Particle Data Group (PDG) 
values [U. The B° -> D^~l + v decays are also used to 
study the systematic uncertainty arising due to the track, 
7T°, and K\ rejections as well as to calibrate the signal 
efficiency. The aforementioned systematic uncertainties 
are estimated by comparing the efficiency before and af- 
ter the application of those vetoes on data and MC. The 
data-MC efficiency ratios for track, 7r°, and K\ vetoes 
are 0.996 ± 0.012, 0.913 ± 0.020, and 1.096 ± 0.020, re- 
spectively. The central values are used to correct the MC 
efficiencies, while the statistical error is treated as a con- 
tribution to the systematic uncertainty. Since the central 
value of the track veto inefficiency is small, no scaling fac- 
tor is applied on the veto efficiency. Instead, the sum of 
the inefficiency and the statistical error is quoted as a 
systematic uncertainty. 




ECL\ 



FIG. 4: .Eecl distribution for doubly tagged events, in which 
one B is fully reconstructed and the other B is reconstructed 
as B° —¥ D^*^£ + u. Points with error bars are data, and the 
solid histogram is a signal MC simulation. 



TABLE I: Summary of the fit result for B° -> D^'l+v 
samples (branching fractions in units of 1CP 3 ). The second 
and third columns show the products of branching fractions, 
where the error on the second column is statistical only. 



Mode 


Measured result 


PDG valucfll] 


B u 


D* 




1.41 ±0.20 


1.34 ± 0.06 


B° 


D* 


~ e + v 


1.62 ±0.18 


1.34 ±0.06 


B° 


D~ 


(Ktttt)h + u 


1.99 ±0.21 


1.98 ±0.12 


B° -> 


D~ 


(Knn)e + v 


1.93 ±0.14 


1.98 ±0.12 


B° -> 


D~ 




0.19 ±0.06 


0.22 ±0.02 


B° -> 


D~ 




0.21 ±0.05 


0.22 ±0.02 



Table Inl lists the signal and background yields for invis- 
ible B decays from the fit while Fig. [5] shows the -Eecl 
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and cos 6 b distributions superimposed with the fit re- 
sult. No significant signal is observed. The signal effi- 
ciency, determined with MC simulations and later cali- 
brated using the doubly tagged B° — > l + v sample, 
is (2.2 ±0.2) x 10 -4 , where the error is dominated by the 
systematic uncertainty. 



TABLE II: Summary of fit yields for the signal and back- 
ground. The normalization of the rare B background contri- 
bution is fixed in the fit. 



Component 


Yield 


Signal 




Generic B background 


131.6t&! 


Non-B background 


9 o 9+21.6 
ZJ - Z -17.0 


Rare B background 


3.7 


Observed events 


121 



The systematic uncertainty associated with the signal 
efficiency is dominated by the -Btag reconstruction effi- 
ciency. The uncertainty on -Btag reconstruction is esti- 
mated by comparing the yield difference between data 
and the corresponding MC sample, generated with a 
proper mixture of generic B and continuum events. The 
Stag yields are extracted by fitting the M^ c distribu- 
tions, and an uncertainty of 8.3% is assigned. System- 
atic uncertainties arising from the requirement of no ad- 
ditional charged tracks nor ir° and K\ candidates are 
estimated to be 1.6%, 2.0%, and 2.0%, respectively, us- 
ing B° — > D^*^£ + iy decays in data. The uncertainty in 
the number of BB pairs is 1.4%. 



TABLE III: Summary of systematic uncertainties arising from 
PDF modeling and components with fixed normalizations. 



Source 


Events 


Signal PDF 


Negligible 


Generic B PDF 


+1.6/- 


1.4 


Rare B PDF 


±0.1 


Rare B fixed yield 


+0.2/- 


0.1 


Non-B PDF 


+1.9/- 


1.3 


Binning effect 


+1.7/- 


1.8 


Sum 


+3.0/- 


2.6 



The uncertainties in the signal yield extraction are 
summarized in Table IIIII The uncertainty due to fixing 
the normalization of the rare B component is obtained 
by varying the rare B yield by the estimated uncertainty 
(±1.9 events). The corresponding variation in the signal 
yield, io'.is ^ s assigned as the systematic uncertainty. For 
each -Eecl PDF, we successively vary the content of each 
histogram bin by ±lcr to obtain a new PDF. The varia- 
tion in the signal yield using the new PDF is calculated 
by performing an unbinned likelihood fit; the quadratic 
sum of all the variations gives the systematic uncertainty 




0.2 0.4 0.6 0.8 1 1.2 

E ECL (GeV) 




COS0 



B 



FIG. 5: The Becl (top) and cos 9b (bottom) distributions 
with fit results superimposed. Points with error bars are data. 
The red cross-hatched region is the signal component on the 
top of the total background shown in the yellow filled his- 
togram. The blue dashed curve is the generic B contribution, 
which is larger than the total because of the negative fit re- 
sult for the non-B background shown in the green dotted 
histogram. The purple hatched area corresponds to the rare 
B contribution. 



for the PDF. The systematic uncertainty arising from 
cos 8b PDFs is negligible. Moreover, the effect of bin 
size is also investigated by choosing different bin sizes 
to model the PDFs. Again, the variation in the signal 
yield is considered as a systematic uncertainty. The to- 
tal systematic uncertainty is computed by summing all 
contributions listed in Table 11111 in quadrature. 

Since there is no significant signal observed, an upper 
limit at 90% confidence level (C.L.) is computed using the 
fit likelihood as a function of the branching fraction. The 
branching fraction is obtained from the signal yield from 
the fit, the signal selection efficiency, and the number of 
BB pairs. The likelihood at each branching fraction is 
obtained using Eq. 1 except that the signal yield is fixed 
in the fit. The systematic uncertainty of the measure- 
ment is taken into account by convolving the likelihood 
function with a Gaussian whose width equals the system- 
atic uncertainty (AS), 
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(B-b') 2 

£ sm ea r (B) = J C$f-j^^d1£!. (2) 

The upper limit on the branching fraction is estimated 
by integrating the likelihood function from zero to the 
bound that gives 90% of the total area. We obtain 
B(B invisible) < 1.3 x 1CT 4 at the 90% C.L. The 
expected upper limit, estimated by applying the same 
method on the MC sample, is 1.1 x 10~ 4 . 

In conclusion, we have performed a search for B — > 
invisible decay with a fully reconstructed B tag on a data 
sample of 657 x 10 6 BB pairs collected at the T(4S I ) res- 
onance with the Belle detector. No significant signal is 
observed, and we set an upper limit of 1.3 x 10~ 4 at the 



90% confidence level for the branching fraction of invisi- 
ble B decay. The limit obtained for B° — > invisible decay 
is the most stringent constraint to date with a hadronic 
tagging method. 
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